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Abstract 

The neutron-drip transition in the dense matter constituting the interior of neutron stars gen¬ 
erally refers to the appearance of unbound neutrons as the matter density reaches some threshold 
density /Odrip- This transition has been mainly studied under the cold catalyzed matter hypothesis. 
However, this assumption is unrealistic for accreting neutron stars. After examining the physical 
processes that are thought to be allowed in both accreting and nonaccreting neutron stars, suitable 
conditions for the onset of neutron drip are derived and general analytical expressions for the neu¬ 
tron drip density and pressure are obtained. Moreover, we show that the nentron-drip transition 
occurs at lower density and pressnre than those predicted within the mean-nucleus approximation. 
This transition is studied numerically for various initial composition of the ashes from X-ray bnrsts 
and superbursts using microscopic nuclear mass models. 
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I. INTRODUCTION 


Born in catastrophic gravitational core-collapse supernova explosions, neutron stars are 
the densest stars known in the universe [l|. According to our current understanding, a 
neutron star contains qualitatively distinct regions. A thin atmospheric plasma layer of light 
elements (mainly hydrogen and helium) possibly surrounds a Coulomb liquid of electrons 
and ions. Below these liquid surface layers, the matter consists of a solid crust made of a 
crystal lattice of fully ionized atoms. With increasing density, nuclei become progressively 
more neutron rich due to electron captures until neutrons start to drip out of nuclei at some 
threshold density pdrip (we shall use the symbol p to denote the mass-energy density). This 
so called neutron-drip transition marks the boundary between the outer and inner crusts. 
The crust dissolves into an homogeneous liquid mixture at about half the density prevailing 
in heavy atomic nuclei. 


The presence of a neutron liquid in the inner crust is expected to play a role in various 
observed astrophysical phenomena such as pulsar sudden spin-ups (so called glitches), and 


the thermal relaxation of transiently accreting neutron stars (see, e.g.. Ref. [2| ). Deter¬ 
mining the onset of neutron drip in dense matter is therefore of utmost importance for the 
modelisation of these phenomena. This transition is generally found to occur at the density 
Pdrip — 4.3 X 10^^ g cm“^ in nonaccreting neutron stars (see, e.g.. Refs. 3l-l6|). This density 
can be shifted due to the accretion of matter from a companion star. For instance, the neu¬ 
tron drip density was found to be given by pdrip — 6.1 x 10^^ g cm“^ (pdrip — 7.8 x 10^^ g cm“^) 
assuming that the ashes of the X-ray bursts consist of pure ®®Fe (^°®Pd respectively) 7, 

On the other hand, these calculations lead to a discontinuous change of the unbound neutron 
density hence also of the neutron chemical potential at the transition. This means that the 
crust is not locally in “chemical” equilibrium: free neutrons will be subject to “chemical” 
forces and will thereby diffuse until the equilibrium is reached. 


In this paper, the condition for the onset of neutron drip is examined considering the 
instability of dense matter against electron captures and neutron emission. Our model of 
dense matter is briefly presented in Section [TTl The usual determination of the neutron-drip 
transition in nonaccreting neutron stars is discussed in Section IIIIl and some misconcep¬ 
tions are pointed out. In particular, we show that the equilibrium nuclei at the onset of 
neutron drip are stable against neutron emission. The stability of dense matter is studied 
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in Section HVl and applications to neutron stars are discussed in Section The cases of 
nonaccreting and accreting neutron stars are considered separately. 


II. MODEL OF DENSE MATTER 


We consider matter at densities high enough that atoms are fully ionized. We further 
assume that the temperature T is lower than the crystallization temperature Tm so that 
atomic nuclei are arranged in a regular crystal lattice. For simplicity, we consider crystalline 
structures made of only one type of ions with proton number Z and mass number A. 
In this case, is given by (see, e.g.. Ref. [l|) 

Tm = ’ ( 1 ) 

e being the elementary electric charge, = (3/(47me))^/^ the electron-sphere radius, rie 
the electron number density, k-Q Boltzmann’s constant, and F^ — 175 the Coulomb coupling 
parameter at melting. Since is generally much lower than the electron Fermi temperature 
dehned by 

^ He — rrieC^ , . 

Ty = , , ( 2 ) 




where fie is the free electron chemical potential, and rUe the electron mass, electrons are 
highly degenerate. Charge screening effects are quite negligible so that electrons will be 
assumed to be uniformly distributed. We will also ignore the small contribution due to 
electron exchange and correlation to the energy density and pressure. The expression of the 
electron energy density Se and pressure Pg in the limit T <C Tp can be found in Chap. 2 of 
Ref. [l|. The free electron chemical potential is given by 






( 3 ) 


where Xr = Ae(37r^ne)^/^ is a dimensionless relativity parameter and Ae = h/{mec) is the 
electron Compton wavelength. For point-like ions embedded in a uniform electron gas with 
number density Ug, the lattice pressure is simply given by 


Pr = 


3 ’ 


where the lattice energy density Sl is given by (see e.g. Chap. 2 in Ref. [l|) 

= Ce^nfZ‘^1^ , 


( 4 ) 


( 5 ) 
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and the crystal structure constant is approximately given by C ~ —1.44 {q] (we have omitted 
here the small contribution due to quantum zero-point motion of ions about their equilibrium 
position). In the ultrarelativistic regime S> 1, the free electron chemical potential and 
the total pressure P = + Pl can be approximately expressed as 


/Xe ~ hc(37r^ne)^/^, 


p 




1 + 


ACaZ^\ 

127r^(hc)^ ’ 

where a = ITlc is the hne structure constant. 


( 6 ) 

( 7 ) 


III. MISCONCEPTIONS ABOUT THE NEUTRON DRIP TRANSITION IN NEU¬ 
TRON STARS 

According to the cold catalyzed matter hypothesis [lO, 11| that will be further discussed 
in Section |Vl the interior of a mature neutron star is in a state of full thermodynamic equi¬ 
librium with respect to all kinds of nuclear and electroweak processes at zero temperature. 
The properties of any layer of the outer crust at some pressure P can thus be determined 
by minimizing the Gibbs free energy per nucleon g, dehned by 

S + P 


9 = 


n 


( 8 ) 


where n is the average nucleon number density, and £ is the average energy density. In 
order to determine the onset of neutron drip, we have to express the Gibbs free energy per 
nucleon for nuclei coexisting with free neutrons and electrons. In the shallowest region of 
the inner crust which we are interested in, the interactions between nuclei and free neutrons 
can be neglected, as shown in Refs. [l2|, ll3| . The average energy density is thus given by 


£ = nxM{A, Z)c^ + £f, + £l + £n — fierrieC^ , 


( 9 ) 


where nx is the number densities of nuclei M{A,Z) their mass (including the rest 
mass of A nucleons and Z electrons) and £n the average energy density of the neutron gas. 
The reason for including the electron rest mass in M{A,Z) is that experimental atomic 
masses are generally tabulated rather than nuclear masses. The last term in Eq.(l9]) is 
introduced to avoid double counting. As pointed in Ref. Q, the concept of ’’free” neutrons 
is ambiguous. In the following we shall consider that a neutron is free if it is unbound 
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in the quantum mechanical sense: its single particle energy in a self-consistent mean-held 
treatment lies above the maximum value of the mean potential. Using the thermodynamic 
identities £e + Pe = '^ehe and £n + Pn = {Pn is the neutron contribution to the pressure, 
Tin the neutron density and /i„ is the neutron chemical potential), the Gibbs free energy per 
nucleon can be written as 


'^X A r7\ 2 I 

g = — Z)cr H- 

n n 


he 


rrieC H- 

3n, 


, Un 

H-hn • 

n 


Let us note that the average nucleon number density is given by 


( 10 ) 


n - A7tx + rin ■ 


(11) 


Moreover, the electric charge neutrality yields 

Ue = Znx ■ (12) 


The Gibbs free energy per nucleon can hnally be expressed as 

4^ 

3 He 

where yn = Un/n denotes the free neutron fraction. In any layer of the outer crust, we 
have |/„ = 0 and P = + Pl- Neutrons start to drip out of nuclei at some pressure 

Tdrip = Pe(^e^‘^) + ^^(ng^'P, Z) whenever 


T Vnl^n y ( 13 ) 


g{A,Z,ne,yn) = (1 - Vn) 


M{A, zy 

A 


+ ^(1 ~ 2/n) 


/ig - mgC 


^(y4, Z, Pdripy dy^) ^y Pdripy Vn 0) 5 


(14) 


or equivalently 


dg 


dyn 




= 0 , 


(15) 


where the partial derivative is evaluated for y^ = 0. This latter condition leads to 


= 0) = j(/l, Z, nf !/„ = 0) = | 


/ie(nf'P) - rUeC^ 


4g^(nfP,Z) 

3 


(16) 

At the neutron drip threshold, we have /i„ ~ (we neglect the small correction of the 

order of a few tens of keV due to neutron-band structure effects [l2, 13|). Equation flTHl) 
can thus be expressed as 


+ -Ce^nyyPzP^ = /if p , 

o 


( 17 ) 
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where 


Z) ^ -M(A,Z)^ + Am„c‘ ^ 


(18) 


Let us note that the equilibrium nucleus at P = Paup must be such as to minimize 
the Gibbs free energy per nucleon. In other words, we must have 


g{A, Z, Pdrip, yn = 0) < g{A', Z\ Pdnp, |/n = 0), 


(19) 


for any values of A' and Z'. In particular, setting Z' = Z and A' = A —AN using Eqs. flT^ . 
ffT6|) . (ITT)) and flTSD . we obtain 

( 20 ) 

emission 


M{A, Z) - M{A - AN, Z) < ANrrin • 
This shows that the nucleus is actually stable against neutron 




A-AN 


X + ANn , 


( 21 ) 


as previously noticed in Ref. 


15|. As can be seen in Table HI the equilibrium nuclei at 


P = -Pdrip as predicted by various microscopic nuclear mass models are indeed stable against 
neutron emission even though they lie beyond the neutron-drip line in the chart of nuclides. 
Let us recall that this line is generally dehned at each value of the proton number Z by the 
lightest isotope for which the neutron separation energy dehned by 


Sn{A, Z) = M{A - 1, Z)c^ - M{A, Z)c^ + , 


( 22 ) 


is negative, i.e., the lightest isotope that is unstable with respect to neutron emission. Be¬ 
cause of pairing and shell effects, many nuclei beyond the neutron-rich side of the neutron 
drip line are actually stable. 

The fact that the equilibrium nucleus at the bottom of the outer crust is stable against 
neutron emission (|2Tll has been often overlooked. The reason may be traced back to the use 
of semi-empirical mass formulae in the early studies of neutron-star crusts. For instance, in 
their seminal work Harrison and Wheeler [w, 11| (the same treatment can also be found 


in the standard textbook from Shapiro and Teukolsky [16|) considered that at high enough 
density the nucleus becomes unstable against neutron emission fl2Tll with AN = 1, at 
which point each nucleus "^X is in equilibrium with its isotope ^~^X and one free neutron, 
as embedded in Eq.(297) of Ref. jll|. In other words, they considered that at the neutron 
drip transition 

M{A,Z)-M{A-l,Z)>mn. (23) 
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TABLE I. Equilibrium nucleus at the bottom of the outer crust of nonaccreting neutron stars, 
as predicted by different microscopic nuclear mass models The corresponding isotope at the 
neutron-drip line in the nuclear chart is indicated. The numbers in parenthesis are the neutron 
separation energies in MeV. The last two columns contain the neutron-drip density and pressure. 


outer crust 

drip line 

Pdrip (g cm“3) 

Rdrip (dyn cm 2) 

HEB-19 ^26Sr (0.73) 

121 Sr (-0.62) 

4.40 X loll 

7.91 X 1029 

HFB-20 i26Sr (0.48) 

i2iSr (-0.71) 

4.39 X loll 

7.89 X 1029 

HFB-21 124si. (0.83) 

i2iSr (-0.33) 

4.30 X loll 

7.84 X 1029 


Using a semi-empirical mass formula, they actually replaced this condition by its continuum 
version 

aM(A,z) 


dA 


> rrin 


The equilibrium nucleus Harrison and Wheeler found at the interface between the outer 
and inner crusts is Z ~ 39 and A ~ 122 (note from Table [I] that this is very close to the 
composition obtained from recent microscopic nuclear mass models). It turns out that this 
nucleus does not satisfy the inequality fl2^ . i.e. it is stable against neutron emission (see 
Appendix). This stems from the fact that the violation of the condition fl2^ corresponds 
to the emission of a huge number of neutrons, one neutron for each nucleus present in the 
crustal layer of consideration. This cannot represent the physical processes associated with 
the neutron-drip transition in neutron star crusts since this would lead to an unphysical 
discontinuous change of the neutron chemical potential. On the contrary, this transition 
can be viewed as the appearance of one free neutron in the whole crustal layer so that the 
neutron fraction tends asymptotically to zero in the thermodynamic limit, and the neutron 
chemical potential thereby varies continuously. It is as if each nucleus were unstable against 
the emission of an inhnitesimally small fraction e of a neutron 

A V . A—e 


X + en. 


Equation 


should thus be replaced by 

M{A, Z) — M(A — e, Z) > errin , 
in the asymptotic limit e ^ 0. It can 


which coincides with Eq. 
continuum approximation leads to Eqs. flTTl) and ffTSjl (see, e.g.. Ref. jl6|). 


(25) 


(26) 


De shown that this 
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Although the previous considerations shed some light on the neutron-drip transition in 
neutron star crusts, they still leave open the question of the actual physical processes that 
can be collectively summarized by Eq. fl25p . This question is of utmost importance for 
accreting neutron stars, whose crusts are generally not in full thermodynamic equilibrium. 


IV. STABILITY OF DENSE MATTER AGAINST ELECTRON CAPTURE AND 
NEUTRON EMISSION 


A. Mean-nucleus approximation 


With increasing density, matter becomes progressively more neutron rich due to the 
capture of electrons, whereby the nucleus transforms into a nucleus z-az'^ with proton 
number Z — AZ and mass number A with the emission of AZ electron neutrino z/g : 


|A + AZe- ^z-Az Y + AZz/g . 


(27) 


The nucleus will be stable against such a process at some given pressure P provided 
the corresponding Gibbs free energy per nucleon is lower than that of the daughter nucleus 
z-Az'^- The expression of the Gibbs free energies per nucleon before and after the electron 
capture can be obtained from Eq. flTOl) . After the electron capture, the average baryon 
number density will be given by 

n"*" = Auy , (28) 


where ny denotes the density of nuclei z-az^- Moreover, the electric charge neutrality 
requires 

7 — A7 

= {Z — AZ)nY = --- vA . (29) 

The Gibbs free energy per nucleon after the electron capture can thus be hnally expressed 
as 


g{A, Z, AZ, n+) 


M{A,Z 


A 


AZ)A ^ Z 


AZ 


A 


, 4£i(n+,Z-AZ) 
he H -;- mgC 

^ 3 nt 


,(30) 


Similarly, the Gibbs free energy per nucleon before the electron capture is simply given by 
g(A, Z, AZ = 0,n~), with 

n~ = Znx = -jn~ , (31) 
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where n~ denotes the average baryon density before the capture. The electron densities n~ 
and are not exactly the same because the pressure P has to remain constant during the 
process. Before the capture, the pressure can be written as 

P = Pe{n-) + PL{n-, Z). ( 32 ) 


After the capture, the pressure can be expressed as 


F = Fe«) + PL(n+,Z-AZ) 


(33) 


The lattice pressure is very small, Pl Pe, so that ^ n+. Let us write = rie + Sue, 
where Ue = n~. Solving Eqs. fl3^ and fl3^ to hrst order in Sue yields 

-1 


6ne = [Piine, Z) - Piiue, Z - AZ)] 


dUf, 


< 0 . 


(34) 


Considering ultrarelativistic electrons and using Eqs. ([6]) and (171), we hnd to first order in a 

Ca 


ISrie 


^2/3 _ 


2/3 


Up <^np. 


(36) 


(3ir2)i/3 _ 

Whereas the electron density varies almost continuously n~ ~ n^, the electron capture will 
be accompanied by a discontinuous change of the baryon density given by 

AZ 


n"*" — n 


( 36 ) 


n Z — AZ 

Expanding the stability condition g{A, Z, AZ = 0,n~) < g{A, Z, AZ,n'^) to hrst order 
in a using Eq. fl5T|) leads to 


Pe + Ce^nl^^ 


Z5/3 - (Z - AZ)5/3 1 3 

AZ 3 




where 




(37) 


(38) 


With further compression of matter, the nucleus may become unstable against the 
capture of electrons accompanied by the emission of free neutrons. Let us consider the 
general case whereby the nucleus transforms into a nucleus z-az^ with proton number 
Z — AZ and mass number A — AN by capturing AZ electrons with the emission of AN 
neutrons n and AZ electron neutrino Ve '■ 


+ AZe- ^ 


A-AN 


Y + ANn + AZz/e 


(39) 
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This process can occur at some pressure P if it leads to a lower value for the Gibbs free 
energy per nucleon. After the electron capture, the average baryon number density is given 
by 

= (A — AN)nY + rin , (40) 

where ny is the density of nuclei z-az'^- Note that AN free neutrons are associated with 
each nucleus so that 

Un = ANny . (41) 

Therefore, the average baryon number density reduces to 

= Any . (42) 

The electric charge neutrality leads to 

n+ = {Z - AZ)ny . (43) 


It follows from these equations that the neutron density can be equivalently expressed as 


n 


n 


AN , 

- 71 

Z-AZ " ’ 


(44) 


for AZ < Z and 




(45) 


for AZ = Z (in this case we must obviously have AN = A and = 0). The Gibbs free 
energy per nucleon after the electron capture and the neutron emission can thus be written 
as 


g{AZ,AZ,AN,nt) = 


Mi A - AN, Z - AZ)c2 AN 

Z-AZ 

A 


A 

, 4S^(nt,Z-AZ) 

fie + X-T- meC 

3 nt 


(46) 


The electron density can be obtained from the requirement that the process occurs at a 
constant pressure P: 


P ^ a(y) + PLin^y 2 - AZ) + F„(n„). 


(47) 


Setting = rig + (frig with Ug = in Eq. fITT)) . and considering AZ < Z, we hnd to hrst 
order in ^rZg 


5ng = [Piirie, Z) - PL{ne,Z - AZ) - P„] 


dPe 

dUf. 


+ 


AN dP„ 
Z — AZ dn„ 


-1 


(48) 
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where dPe/due has to be evaluated at the density Ug, whereas Pn and dPn/drin at the density 
ngAN/(Z—AZ). Typically dPn/dun -C dPe/dng. We shall therefore neglect the derivative of 
the neutron pressure in Eq. fl48p . Expanding the stability condition g{A, Z, AZ = 0, AN = 
0, n~) < g{A, Z, AZ, AN, to hrst order in a using Eqs. flTB|l and (1451) . we hnd 


/ie + Ce^nl^^ 
AN 


ZV3 _ _ AZfA ^ 


AZ 


{l^n 


rrinC 


AZ 





(49) 


where 




M{A - AN, Z - AZy - M{A, Zy + rUnC^AN 

AZ 


+ mgC^, 


(50) 


and we have assumed AZ > 0. We have also neglected the term arising from the expansion 
of the neutron chemical potential since 


AN dgin/ 

Z — AZ d/ig/dug 


dPn/dUn 

dPg/dng 


(51) 


In case of neutron emission without any electron capture (i.e. AZ 
hnd 


M(A, Z) - Mi A - AN, Z) < AN^ . 


0 and AN > 0), we 
(62) 


It is to be understood that pn and its derivative are evaluated at density ngAiV/(Z — AZ). 
Therefore, fin > whenever AN > 0. In the limiting case AZ = Z, the stability 

condition takes a form similar to Eq. fll7p except that in Eq. ffTSD is now replaced by 

/X„(n+) > TTlnC^. 

This analysis thus leads to a discontinuous change of the neutron density hence also of 
the neutron chemical potential (unless of course AN = 0). Moreover, as can be seen by 
comparing Eq. fl52l) to Eq. fl20p . the nuclei are not necessarily stable against neutron 
emission contrary to the general considerations of Section HTTl These unphysical results stem 
mainly from our implicit assumption (arising from the mean-nucleus approximation) that 
all nuclei in a given crustal layer (i.e. at a given pressure P) will emit neutrons. In reality, 
only some fraction of nuclei may become unstable against neutron emission so that the crust 
will be most presumably composed of an admixture of nuclei and 


11 










B. Beyond the mean-nucleus approximation 


Let us consider that after the onset of electron capture and neutron emission pro¬ 
cesses fl5^ . only some fraction of nuclei will transform into case, the 

Gibbs free energy per nucleon can be written as 


g = —MiA, Z)c^ + —MiA - AN, Z - AZy + — 
n n n 


Tfl^C 


n. 


+ fi'L H- fj-n 5 (53) 


n 


where ql denotes the lattice contribution, for which we use for simplicity the linear mixing 
rule (see, e.g.. Section 2.4.7 in Ref. [fl). The lattice energy density is thus given by 


S 


m 

L 


Z-yAZ 


il-y)Z^A + y{Z-AZf/^ , 


(54) 


where y = ny/{nx + ny). For the lattice contribution to the pressure, we hnd P™ = T™/3. 
The lattice term can hnally be expressed as 


9l 


+ Pl 


n 


4^ 

3 n 


(55) 


The baryon number density is given by 


n = Anx + [A — AN)ny + . 

Free neutrons are associated with nuclei therefore 

Un = ANny . 


(56) 


(57) 


The baryon density can thus be expressed as 

n = A{nx + ny). (58) 

The electric charge neutrality requires 


He = Znx + {Z — AZ)ny . 


(69) 


Using these equations, the electron and neutron densities can be equivalently written as 

(60) 

(61) 


Z-yAZ 
fie = - J - n . 

yAN yAN 
Un = - —n = — - —He . 


A 


Z-yAZ 
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Finally, we obtain for the Gibbs free energy per nucleon 


g{A, Z, AN,AZ, rie.y) 


, ,M{A,Zy M(A- AN,Z - AZy 

( 1 -y) ^ . +y- ^ ^ 


A 


A 


Z-yAZ 

A 

-^ceyy 

3A 


/^e Tfl^C 


+ 


AN 


A 


-yy-r 


{l-y)z^/^ + y{Z - Azy^ 


( 62 ) 


Note that y can be treated as a continuous variable, which measures the degree of neutron- 
ization of matter. Before the onset of electron capture and neutron emission, y = 0 so that 
= 0 and 

n- = ^n- . (63) 


The pressure is given by Eq. fl32|) . and the Gibbs free energy per nucleon is given by 
g{A, Z, AN^ AZ, ?/ = 0). With increasing pressure, some inhnitesimally small fraction dy 
of nuclei may become unstable against neutron emission leading to an infinitesimally small 
neutron density 


dUn 



( 64 ) 


This shows that the neutron density varies continuously across the transition contrary to the 
mean-nucleus approach discussed in the previous section. For the nucleus to be stable 
against electron captures and neutron emissions, the following condition must be satished: 


g{A z, AN, AZ, n+, dy) > g(A, Z, AN, AZ, 7i^,y = 0), 


(65) 


where the electron densities and are related by 


P = Pfin;) + PLiy, Z) = Pfint) + pyy, Z, AZ, dy) + P„(dn„). 


( 66 ) 


This equation is similar to Eq. fl47)l obtained in the mean-nucleus approximation, except that 
the neutron pressure term is now vanishingly small and can thus be dropped. Indeed, using 
Eq. fl6T)) we have 

Pnidun) = Pn , (67) 

where dPn/dun is evaluated at = 0. On the other hand, the many-body theory for a 


dilute neutron gas with scattering length a < 0 yields 


17| 


8n o h^7r|a|n„ 

— mnC -- 


Ur. 


5 2mr, 


mr. 


1 - 


357r 


(11 - 21n2)A:F|a| 


( 68 ) 
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so that dPri/dun —)■ 0 in the limit = 0. Therefore, Pn{dnn) = 0. 

Setting = He + Srie with Ug = n~ in Eq. (]66|l . we hnd to hrst order in Sug 

-1 


5ng = [PL{ng, Z) - Z, AZ, dy)] 


dPg 

dug 


( 69 ) 


Expanding Eq. fl^ to hrst order in a using Eq. (1^ and (1^ . we hnally obtain after some 
algebra 


/ie + Ce^nl^^ 


Z5/3 _ 1 3 

AZ 3 




(70) 


where 




_ M{A - AN, Z - AZy - M{A, Zy + my AN 


AZ 


+ mgm , 


(71) 


and we have assumed AZ >0. In case of neutron emission without any electron capture 
(i.e. AZ = 0 and AN > 0), we hnd 


M{A, Z) - M{A - AN, Z) < ANm^ • 


(72) 


Note that this last inequality is independent of the electron background, and coincides with 
the stability condition (120|) . In other words, a nucleus unstable against neutron emission in 
vacuum will be also unstable in dense matter, and therefore such a nucleus cannot exist in 
equilibrium. This means that the neutron-drip transition in dense matter must be triggered 
by electron captures. 


V. MATTER NEUTRONIZATION IN NEUTRON STARS 


A. Nonaccreting nentron stars 


We consider here solitary neutron stars formed in supernova explosions of single massive 
stars (with a mass M > 8Mq, Mq being the mass of the Sun). Initially, the newly born 
neutron star is very hot (T ~ 10^*^ K) and fully huid. This is the “hot scenario” of crust 
formation. A few days after the neutron star birth the outer layers are still so hot (T > 
2 X 10® K, see e.g. Fig. 34 in Ref. 0 ) that nuclear reaction rates are sufficiently high to 
keep matter close to the nuclear equilibrium corresponding to the minimum of the Gibbs free 
energy per nucleon g at given temperature T (decreasing) and pressure P (increasing). The 
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star cools by neutrino emission and after a few months its outer layer crystallizes forming a 
solid crust beneath an ocean of a fluid hot plasma. Following Refs. [lO|,lllt], it is assumed that 
during further cooling the crust remains in nuclear equilibrium until it eventually becomes 
cold and fully “catalyzed”. In order to reach this state, all possible kinds of electroweak 
and nuclear reactions should have sufficient time to be completed. This means that electron 
captures and neutron emission processes fl3^ with all possible values of AZ and AN should 
be considered. The neutron drip transition can thus be found by minimizing Eq. fl50|) with 
respect to AZ and AN . Obviously, the number AZ of electrons that can be captured by 
a nucleus must be lower than Z, and the number AN of neutrons that can be emitted 
must be lower than A. Moreover, the number of bound neutrons in the daughter nucleus 
z-Az^ should be positive, thus leading to AZ > Z + AN — A. On the other hand, we have 
previously shown that the neutron emission must be accompanied by electron captures so 
that AZ > 0 therefore AN > A — Z. All in all, we have 


Z>AZ>Z + AA-A, A> AN > A - Z . 


(73) 


Using Eqs. ([6]), fl62ll . (1701) . and flTTl) . the Gibbs free energy per nucleon at the neutron drip 
transition is given to lowest order in a by 


g{A,Z,AN,AZ,ne,y = 0) = 


M{A, zy 


+ 


A 

A(37r2)V3 _ 

ZAN 2 

rrinC 


AZ 


+ 


Z5/3 ( 1 _ 


AZ 


Z M{A- AN,Z - AZy 
AZ A 

Z(Z-AZ)^/3 


+ 


AZ 


AAZ 


(74) 


The lowest value of the Gibbs free energy per nucleon, g = my, is thus obtained for 
AZ = Z and AN = A: the electron captures proceed until the complete disintegration of 
the nuclei. Similar considerations were previously discussed in the seminal work of Ref. j^. 
It can be easily seen that in this case Eq. (|7nil reduces to Eq. (ITTll . Gonsidering that electrons 
are ultrarelativistic using Eqs. (jb]) and ([7]), the baryon density and pressure at neutron drip 
are approximately given by 

^drip(A, Z) —^7 n/* -.J 1 + 


-Pdrip(A, Z) 

with Z) given by Eq. 


A /if P(A, Z)3 

z ?,yncf 

/if P(A, Z)4 


ACa 


1 + 


-Z^A 

(8l7r2)V3 

ACaZ"^/^^-^ 


-3 


127r2(hc)3 (8l7r2)^/3 


(75) 

(76) 
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B. Accreting neutron stars 


The accretion of matter onto a neutron star from a stellar companion may change the 
initial constitution of its crust, which was originally formed in the “hot scenario” described 
in the previous section. For an accretion rate M = 10“® MQ/yr the original outer crust 
is replaced by accreted matter in 10^ yr, while the accretion stage in the low-mass binary 
systems can last for 10® yr. At densities above ~ 10® g cm“® matter is highly degenerate 
and relatively cold (T ^ 5 x 10® K) so that thermonuclear processes are strongly suppressed; 
their rates are many orders of magnitude lower than the compression rate due to accretion. 
This is the “cold scenario” of crust formation. On the other hand, the matter composition 
can still be altered due to electron captures, neutron emissions, and at high enough densities 
pycnonuclear reactions. Multiple electron captures are unlikely and need not be considered. 
For example, the double electron capture 2e“ -|- ®®Fe —>■ ®®Mn -f- 2ve occurs on a timescale 
of 10^® yr 19|. The onset of neutron drip can thus be determined from Eqs. (I70|) and flTTj) 
with AZ = 1 and AA > 0. In the limit of ultra relativistic electrons, using Eq. (|6]), the 


average baryon density for the onset of electron captures is approximately given by 


^drip-acc (A,Z) 


A/if 


1 + 


Ca 


_i^Z- 1)®/® + 


(77) 


Z 37r®(hc)® [ ' 

^drip-acc^^^ Z) = M{A -AN,Z- l)c^ - M{A, zy + ANrrinC^ + rrieC^ . (78) 

The threshold pressure Pdnp-acc can be obtained from Eq. ([7]) with /ie obtained from Eq. (170]) : 


P. 


drip-acc 


(A,Z) 


^drip-accZ)4 

127r2(hc)® 


1 + 




(8l7r2)V3 


1 + 


Ca 


(37r2)i/3 


^ 5/3 - 1 )®/® + 


^2/3 


-4 


(79) 


For the sake of comparison, let us calculate the neutron drip transition in accreting 
neutron stars using a strict mean-nucleus approximation. For this purpose, we will apply 
the Mackie-Baym compressible liquid drop model widely used in previous calculations of 
accreted neutron-star crusts [ 7 , 8, 201. Assuming that the ashes of the X-ray bursts consist 
of pure ®®Fe, the equilibrium nuclei were found to become progressively more neutron rich 
with increasing density due to electron captures. At density 5.65 x g cm“® corresponding 
to the pressure 1.23 x 10®® dyn cm“^, the nucleus ®®Ar is found to be unstable against one 
neutron emission, 

®®Ar + e" Cl + n + z/g . (80) 
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Solving numerically Eq. flTOj) with AN = 1 and the masses calculated with the same 
compressible liquid drop model as in Refs. [?, 8| in the absence of a neutron gas, we 
now hnd that the neutron drip transition actually occurs at a lower density, given by 
Pdrip-acc — 5.09 X10^^ g cm“^. The neutron-drip pressure is now given by Pdrip-acc — 1-07x 10^° 
dyn cm“^. For ashes made of ^°®Pd, using the Mackie-Baym compressible liquid drop model 
of accreted crust, the neutron drip nucleus ^°®Ge is found to decay into ^°^Ga, 


106 / 


Ge -I- e —Ga + 3n + Ue , 


(81) 


at density 6.66 x 10^^ g cm“^ (pressure 1.38 x 10^° dyn cm“^). Solving now numerically 
Eq. fl70]l with AN = 3, the neutron-drip density decreases to Pdrip-acc — 5.31 x 10^^ g cm“^. 
The drip pressure is Pdrip-acc — 1-02 x 10^° dyn cm“^. Let us remark that the nuclei ^®G1 and 
^°^Ga are unstable and undergo further electron captures accompanied by neutron emissions. 
We found that the errors of the analytical expressions (1771) and (1791) lie below about 0.1%. 
Approximating the neutron drip density by Pdrip-acc ~ J^’T-drip-acc where m is the unihed 
atomic mass unit, and using Eq. (ITTj) . leads to an error of about 0.7%. 

It should be remarked that our values of the neutron-drip density calculated within the 
mean-nucleus approximation differ from those previously obtained in Refs. M, Isj because of a 
different treatment of neutrons produced by electron captures. In Refs. Qff these neut/ons 
were kept inside the nucleus, assuming that free neutrons outside remained unaffected. 
The daughter nucleus was then found to be unstable against neutron emission alone. In 
the present approach, we minimize the Gibbs free energy per nucleon (for a given atomic 
number) without any further constraint. This new procedure yields values for the neutron- 
drip density and pressure that are closer to the exact results obtained from Eq. (ITOj) . Indeed, 
using the new version of the Mackie and Baym model of accreted crust for the ®®Fe ashes, 
we are now getting 5.7 x 10^^ g cm“^ instead of 6.1 x 10^^ g cm“^ obtained in Refs. { 7 , 8 ]. 


In the case of the ^°®Pd ashes the di 


7.8 X 10^^ g cm ^ obtained in Refs 


.0 


erence is even larger, 6.7 x 10^^ g cm ^ instead of 

8 |. 


As the nuclei from the ashes of X-ray bursts sink into the crust, their proton number Z 
decreases due to electron captures whereas A remains unchanged (we assume that pycnonu- 
clear reactions may only occur in the inner crust, where nuclei are immersed in a neutron 
liquid). At some point, the daughter nuclei will be so neutron rich that neutrons will be 
emitted. Therefore, the neutron-drip transition will occur when the threshold electron chem- 
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ical potential /if"" for neutron emission (i.e. AN > 0 and AZ = 1) will become lower than 
the threshold electron chemical potential /rf for electron capture alone. Using Eqs. fl38|l and 
the condition Z) < Z) yields 


S^Nn{A, Z-l)= M{A - AN, Z-1)- M{A, Z - 1) + AiVm„ < 0 . 


(82) 


In other words, the nucleus marking the neutron drip point is such that the nucleus 
z-iX is unstable against neutron emission. Consequently, for any given value of the mass 
number A, the proton number Z of the nuclei present in the outer crust will decrease until 
the AiV-neutron separation energy SANn{A,Z — 1) becomes negative: this will mark the 
onset of neutron drip. 

We have determined in this way the neutron-dripping nucleus in the crust of accreting 
neutron stars. As for the initial composition of the ashes, we considered two different 


scenarios: hrst, the ashes are produced by an rp-process during an X-ray burst 


second, the ashes are produced by steady state hydrogen and helium burning [22| as expected 


21 |. and 


to occur during superbursts 


23l |. We have calculated the neutron-drip density and pressure 
by solving numerically Eq. fTfOl) without any further approximation considering all possible 
neutron emission processes. Results are shown in Tables H H IIIII and IIVI for the three different 
microscopic nuclear mass models HFB-19, HFB-20 and HFB-21, which are based on the self- 
consistent Hartree-Fock-Bogoliubov method 2^. The neutron drip density and pressure in 
accreting neutron stars are thus found to be quite sensitive to the composition of the ashes. 
In reality, the crust of an accreting neutron star is expected to contain an admixtures 
of various nuclides. Therefore, the onset of neutron drip will be determined by the most 
unstable nucleus. Depending on the adopted mass model, the threshold density is thus found 
to vary from 2.80 x 10^^ g cm“^ to 6.13 x 10^^ g cm“^ (Tdrip-acc ranging from 4.72 x 10^® 
dyn cm“^ to 13.1 x 10^® dyn cm“^) for ordinary burst ashes, and from 3.23 x 10^^ g cm“^ 
to 6.13 X 10^^ g cm“^ (Udrip-acc ranging from 6.11 x 10^® dyn cm“^ to 13.1 x 10^® dyn cm“^) 
for superburst ashes. For comparison, as indicated in Table [H the neutron-drip density in 
catalyzed matter lies in the range 4.30 — 4.40 x 10^^ g cm“^ (Tdrip ranging from 7.84 x 10^® 
dyn cm“^ to 7.91 x 10^® dyn cm“^). This analysis shows that the neutron-drip transition in 
accreting neutron stars may occur either at a higher or at a lower density than in nonaccreting 
neutron stars. In tables IVl and IVIl the predictions from the microscopic nuclear mass models 
are compared to those from the Mackie and Baym compressible liquid drop model employed 
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TABLE II. Neutron drip transition in the crust of accreting neutron stars, as predicted by the HFB- 
19 microscopic nuclear mass model: mass and atomic numbers of the dripping nucleus, number 
of emitted neutrons, density and corresponding pressure. The mass numbers A are listed from 
top to bottom considering that the ashes are produced by ordinary X-ray bursts (upper panel) or 
superbursts (lower panel). See text for details. 


A 

Z 

/?drip—acc 

(10^^ g cm ^) 

Rdrip-acc (10^® dyn cm“2) 

104 

32 

1 

4.76 

9.08 

105 

33 

1 

3.31 

5.74 

68 

22 

1 

4.00 

7.78 

64 

20 

3 

5.74 

12.0 

72 

22 

1 

5.15 

10.1 

76 

24 

1 

4.97 

10.0 

98 

32 

1 

3.34 

6.13 

103 

33 

1 

2.82 

4.76 

106 

34 

1 

3.65 

6.73 

66 

22 

1 

3.46 

6.68 

64 

20 

3 

5.74 

12.0 

60 

20 

1 

3.23 

6.11 


in Refs. j3, 
as in Refs. 


8 |. For this purpose, we have considered the same initial composition of ashes 


7 




8 |. For ^®Fe ashes, all models predict the same atomic number Z = 18 for the 
dripping nucleus and the same number AN = 1 of emitted neutrons. On the other hand, 
the HFB models predict substantially lower values for the average neutron drip density and 
pressure. Much larger deviations are found in the case of ^‘^®Pd ashes, mainly due to different 
predictions for the dripping nucleus and the number of emitted neutrons. In particular, the 
neutron-drip densities and pressures, as obtained from HFB-21 (which is favored over the 
two other HFB models for the reasons given in Refs. Hfl) and the compressible liquid 
drop model, differ by almost a factor a two. This suggests that t 
crusts could be substantially different from those found in Refs. 


le composition of accreted 

asn 


,l20|. 


19 












TABLE III. Same as Table HIl for the HFB-20 microscopic nuclear mass model. 


A 

Z 

/?drip—acc 

(10^^ g cm 

-Pdrip-acc (10^® dyn cm“2) 

104 

32 

1 

4.82 

9.23 

105 

33 

1 

3.38 

5.91 

68 

22 

1 

4.09 

8.01 

64 

20 

3 

6.13 

13.1 

72 

22 

1 

5.21 

10.3 

76 

24 

1 

5.01 

10.1 

98 

32 

1 

3.34 

6.14 

103 

33 

1 

2.80 

4.72 

106 

32 

1 

5.26 

10.1 

66 

22 

1 

3.54 

6.88 

64 

20 

3 

6.13 

13.1 

60 

20 

1 

3.27 

6.20 


TABLE IV. Same as Table [E for the HFB-21 microscopic nuclear mass model. 


A Z AN Pdrip-acc (10^^ g cm 3) Pdrip-acc (10^^ dyn cm 2) 


104 

32 

1 

4.85 

9.31 

105 

33 

1 

3.42 

6.01 

68 

22 

1 

4.13 

8.12 

64 

20 

3 

5.84 

12.3 

72 

22 

1 

5.35 

10.6 

76 

24 

1 

5.02 

10.2 

98 

32 

1 

3.42 

6.33 

103 

33 

1 

2.83 

4.79 

106 

34 

1 

3.65 

6.72 

66 

22 

1 

3.58 

6.98 

64 

20 

3 

5.84 

12.3 

60 

20 

1 

3.36 

6.43 
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TABLE V. Neutron drip transition in the crust of accreting neutron stars, as predicted by dif¬ 
ferent nuclear mass models for ^®Fe ashes: atomic number Z of the dripping nucleus, number of 
emitted neutrons, density and corresponding pressure. Three different microscopic Hartree-Fock- 
Bogoliubov nuclear mass models (HFB) are compared to the compressible liquid drop model of 
Mackie and Baym (MB). See text for details. 



HFB-19 HFB-20 HFB-21 

MB 

z 

18 

18 

18 

18 

AN 

1 

1 

1 

1 

/^drip—acc (10 S ) 

4.49 

4.50 

4.38 

5.09 

Thrip-acc (10^® dyn cm-2) 

9.02 

9.06 

8.74 

10.7 

TABLE VI. Same as 

Table El fo 

r ^^®Pd ashes. 



HFB-19 HFB-20 HFB-21 

MB 

Z 

34 

32 

34 

32 

AN 

1 

1 

1 

3 

/^drip—acc (10 g CIXL ) 

3.65 

5.26 

3.65 

5.31 

Thrip-acc (10^® dyn cm-2) 

6.73 

10.1 

6.72 

10.2 


VI. CONCLUSIONS 


With increasing pressure, dense matter in the outer crust of a neutron star becomes 
progressively more neutron rich. At some point, neutrons start to drip out of nuclei. In 
both accreting and nonaccreting neutron stars, the nucleus at the onset of neutron drip is 
shown to be stable against neutron emission, but is actually unstable against electron cap¬ 
ture accompanied by neutron emission. After examining the occurrence of such processes 
in both accreting and nonaccreting neutron star crusts, we have obtained general analytical 
expressions for the density and pressure at the onset of neutron drip considering that only 
a vanishingly small fraction of nuclei become unstable. In this way, we have shown that 
the spurious discontinuous change in the density of unbound neutrons that was found in 
previous studies of accreting neutron stars {7,8] arises from the use of the mean-nucleus ap¬ 
proximation. As a consequence, this approximation overestimates the neutron-drip density 
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and pressure. 

We have also studied numerically the transition between the outer and inner crusts of 
accreting neutron stars for various initial compositions of the ashes. In particular, we con¬ 
sidered ashes produced from both X-ray bursts and superbursts. For this purpose, we have 
made use of experimental atomic mass data complemented with microscopic atomic mass 
tables. The neutron-drip density is found to be shifted to either lower or higher values than 
in nonaccreting neutron stars, depending on whether accreting neutron stars exhibit ordi¬ 
nary bursts or superbursts. The crust of accreting neutron stars is also predicted to contain 
various ultradrip nuclei so that it is necessary to extend the calculations of atomic masses 
beyond the neutron drip line. Finally, large deviations were found between the predictions 
of microscopic mass models and the more phenomenological liquid drop models employed 
in Refs. j3,1^ , thus suggesting that the composition and the properties of accreted neutron 
stars crusts could differ substantially from those predicted in Refs. |7|, l 8 |. 


Appendix A: Neutron drip transition in cold catalyzed matter with the Harrison- 
Wheeler model 


In the seminal work from Wheeler and collaborators jld. Illl] (see also Ref. [l^), nuclear 
masses were calculated using the semi-empirical formula of Green 


M(A, Z) = rriit 


biA + b2A'^/^ - bsZ + b^A { - - j\ 


(Al) 


where bi = 0.992064, 62 = 0.01912, 63 = 0.00084, 64 = 0.10178, 65 = 0.000763, and 
niu = M(^®0)/16. The equilibrium nucleus present in the crust of a nonaccreting neutron 
star can be determined by minimizing the Gibbs free energy per nucleon g with respect to 
both Z and A for any given pressure P. Neglecting the lattice contribution to q and treating 


111 


Z and A as continuous variables lead in particular to the following relation 

The actual values of Z and A can be found from the beta-equilibrium condition [n| 

dM{A, Z) 


dZ 


= nip 


C , 


(A3) 


where the free electron chemical potential is obtained from the requirement that the electron 
pressure be given by P. The neutron-drip transition occurs at some pressure Pdrip such that 
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the following condition is fulfilled: 


dM{A,Z) 


= 


(A4) 


Solving Eq. flA4p using Eqs. flA2p and flAip yields Zdrip ^ 39.09 and A 
adopted the same value for the neutron mass as in Ref. [n|, namely m 
dripping nucleus is stable against neutron emission: 


drip — 122.0. We have 
n = 1.008982m,r. The 


M(Adrip, ^drip) - Af(Adrip - 1, ^drip) - rUn ^ “8.6 X 10 ^rUu < 0 . (A5) 
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